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We perform a detailed phenomenological analysis of how well hadronization in nuclear environ- 
ments can be described in terms of effective fragmentation functions. The medium modified frag- 
mentation functions are assumed to factorize from the partonic scattering cross sections and evolve in 
the hard scale in the same way as the standard or vacuum fragmentation functions. Based on precise 
data on semi- inclusive deep-inelastic scattering off nuclei and hadron production in deuteron-gold 
collisions, we extract sets of effective fragmentation functions for pions and kaons at NLO accuracy. 
The obtained sets provide a rather accurate description of the kinematical dependence of the ana- 
lyzed cross sections and are found to differ significantly from standard fragmentation functions both 
in shape and magnitude. Our results support the notion of factorization and universality in the 
studied nuclear environments, at least in an effective way and within the precision of the available 
data. 
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I. INTRODUCTION AND MOTIVATION 

In spite of the remarkable phenomenological success 
of Quantum Chromodynamics (QCD) as the theory of 
strong interactions, a detailed quantitative understand- 
ing of hadronization processes is still one of the great 
challenges for the theory. Hadronization is the mecha- 
nism by which a final-state quark or gluon, excited in 
some hard partonic interaction, dresses itself and devel- 
ops into an observed final-state hadron. As such, it is 
sensitive to physics happening at long distances and time 
scales, and perturbative tools that have earned QCD its 
present standing, simply fall short. On the other hand, 
studying hadronization process is a path to understand 
the phenomenon of confinement and opens a window to 
the non-perturbative domain of QCD. 

The last few years have seen a significant improve- 
ment in the perturbative QCD (pQCD) description of 
scattering processes with identified produced hadrons. 
More specifically, precise determinations of fragmen- 
tation functions for various different hadron species 
have been performed P, 0, 0j- Fragmentation func- 
tions parametrize the non-perturbative details of the 
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hadronization process and, by virtue of the factorization 
theorem, are assumed to be universal, i.e., process in- 
dependent Q. One of the most interesting results of 
these studies is that the standard pQCD framework not 
only reproduces the data on electron-positron annihi- 
lation into hadrons with remarkable precision, but de- 
scribes equally well other processes, like semi-inclusive 
deep-inelastic scattering (SIDIS) and hadron production 
in proton-proton collisions Using the wealth of data 
with identified hadrons in a global QCD analysis not only 
increases the precision and kinematical coverage of frag- 
mentation functions, but more importantly, supports the 
ideas of factorization and universality in the kinematical 
domain accessed by these experiments. These concepts 
are the starting points for the successful pQCD descrip- 
tion of hard scattering processes. 

For more than thirty years, it is known that results 
for hadron production processes occurring in a nuclear 
environment can differ significantly from similar exper- 
iments involving only light nuclei or proton targets [f|. 
More recently, detailed kinematic distributions for pro- 
duction rates of identified hadrons in SIDIS off different 
nuclei have been provided by the HERMES experiment 
Q. The origin of the observed differences induced by 
the nuclear media, has been attributed to a variety of 
conceivable mechanisms besides the well-known modifi- 
cation of parton densities in nuclei. Ideas range from 
interactions between the nuclear medium and the seed 
partons before the hadronization takes place, as imple- 
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mented through pQCD inspired models, to interactions 
between the medium and the final-state hadrons, usually 
formulated in a nuclear and hadron physics language; for 
recent reviews, see Refs. 0j@|- Most models reproduce, 
with different degree of success, some features of the data, 
in spite of very different, even orthogonal theoretical ap- 
proaches and ingredients. 

In the case of the so called initial-state nuclear effects, 
such as those observed in inclusive deep-inelastic scat- 
tering (DIS) off nuclei and the Drell-Yan process on 
nuclear targets [10|, it has been shown that the influ- 
ence of the nuclear environment can be effectively fac- 
torized into a set of modified nuclear parton distribu- 
tion functions (nPDFs). The nPDFs scale in energy with 
standard evolution equations, and conventional partonic 
hard scattering cross sections can be utilized in calcu- 
lations. Within the precision of the available data, this 
approach has been demonstrated to be an excellent ap- 
proximation at next-to-leading order (NLO) accuracy of 
pQCD, and allows one to continue to exploit the key 
features of factorization and universality in a nuclear en- 
vironment [IH EH EH • From the QCD point of view, the 
nPDFs carry all the non-perturbative information related 
to the probed nuclei, just as the standard parton densities 
parametrize our knowledge of the nucleon. Both need to 
be obtained in global QCD fits to data or approximated 
within a given non-perturbative model. 

It is quite natural then to ask if the idea of factorization 
can be extended to final-state nuclear effects and to ex- 
plore how good such an approximation works in practice 
or, alternatively, to determine where and why it breaks 
down. From a theoretical point of view, however, the an- 
swer to the question whether the introduction of medium 
modified fragmentation functions (nFFs) is a viable con- 
cept is not obvious. On the one hand, interactions with 
the nuclear medium may spoil the requirements of fac- 
torization theorems, but on the other hand any estimate 
of the factorization breaking effects is strongly model de- 
pendent. 

Assuming that factorization holds, the medium modi- 
fied fragmentation functions should contain (factorize) all 
the non-perturbative details relevant for the computation 
of hard processes with identified hadrons and would be 
interchangeable from one process to another, i.e., univer- 
sal. On top of that, such a theoretical framework is well 
defined, has predictive power, and can be systematically 
extended beyond the leading order (LO) approximation. 
The inclusion of higher order QCD corrections and the 
possibility to use different observables simultaneously in 
a global analysis have both been proven to be crucial for a 
successful extraction of vacuum fragmentation functions 
(FFs) All these aspects can be explicitly tested also 
in a nuclear environment, using data from an increasing 
number of experiments that have performed very precise 
measurements of inclusive hadron production off nuclear 
targets. They comprise semi-inclusive measurements of 
hadron multiplicities by HERMES [6|, as well as exper- 
iments in deuteron-gold (dAu) collisions at the BNL- 



RHIC 0,0,13. Both processes show clear signals 
of a non-trivial nuclear dependence in the hadronization 
mechanism. 

In addition to our primary goal, which is testing the 
factorization properties of nFFs in a consistent theoreti- 
cal framework, we also aim to constrain them as precisely 
as possible from the different data sets currently avail- 
able. This will allow one to compare our results for the 
nFFs with the different model estimates and mechanisms 
proposed to model hadronization in a nuclear medium 
0, @| • Together with our analysis of vacuum FFs [l[ , the 
present work also serves as baseline for ongoing studies 
of processes with detected hadrons in heavy ion colli- 
sions performed at BNL-RHIC and at the CERN-LHC 
in the future. These projects are aimed at investigating 
the properties of hot and dense QCD matter, and require 
precise knowledge on hadronization under both normal 
and such extreme conditions 

The paper is organized as follows: in the next Section, 
we very briefly summarize the pQCD framework for FFs 
and discuss how to extend it to nFFs, specifically for 
the processes on which we will base our analysis. We 
also review the different aspects of the data that suggest 
medium induced effects in fragmentation processes. In 
Sec. Ill, we first show how the main features of the exper- 
imental results can be reproduced in a very basic scheme 
of medium modifications, associated with a simple pa- 
rameterization of nFFs which can be motivated by intu- 
itive arguments. Next, we present our results of a more 
refined determination of the medium modified fragmen- 
tation functions for pions at NLO accuracy. We briefly 
comment on different centrality classes in deuteron-gold 
collisions at RHIC. Finally, we extend our analysis to 
measurements of kaon production in a nuclear environ- 
ment. We summarize our results in Sec. IV. 



II. QCD FRAMEWORK FOR MEDIUM 
MODIFIED FRAGMENTATION FUNCTIONS 

A. Basic properties of FFs 

In the naive parton model, fragmentation functions 
Df 1 (z) are simply taken as the probabilities for a final- 
state parton of a given flavor i to produce a specific 
hadron H, carrying a fraction z of its four-momentum. 
This intuitive picture can be consistently extended to 
the field-theoretical language of QCD by defining frag- 
mentation functions in terms of bifocal operators in a 
certain factorization scheme [l8jj. The most common 
choice for the latter is the MS scheme, which we also 
adopt throughout this work. Along with analogous op- 
erator definitions for the parton distribution functions 
(PDFs), the factorization of short- and long-distance con- 
tributions to one-hadron inclusive cross sections is then 
precisely determined, and higher order QCD corrections 
can the be systematically taken into account. While the 
non-perturbative but universal FFs (and PDFs) need to 
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be extracted from data through global QCD analyses, 
both the relevant hard partonic scattering cross sections 
and the scale evolution of FFs (and PDFs) are calculable 
within pQCD. Corrections to this factorized framework 
emerge as inverse powers of the hard scale characterizing 
the one-hadron inclusive cross section, like the virtuality 
Q 2 of the exchanged virtual photon in e + e~ annihila- 
tions or the transverse momentum px of the produced 
hadron in proton-proton collisions. Provided that Q 2 or 
Pt are large enough, these higher twist contributions can 
be safely neglected in phenomcnological applications. 

The scale dependence of the Df is governed by a set 
of coupled renormalization group equations very similar 
to those for PDFs. They schematically read 
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is the matrix of the singlet timelike evolution kernels. 
The symbol <g> denotes a standard convolution. The NLO 
splitting functions Pf?^ have been computed in [111, HO] or 
can be related to the corresponding spacelike kernels by 
proper analytic continuation [2lj | . Recently, the diagonal 



splitting functions Pq q and Pgg ' have been calculated 
up to next-to-next-to- leading order accuracy [22} ■ 

Apart from the requirement of a sufficiently large hard 
scale to suppress power corrections, the applicability of 
FFs is in practice further limited to the range of medium- 
to-large momentum fractions z [3, Hi ■ While large per- 
turbative contributions to the splitting kernels PgP and 

frp\ 

Pgg at small z can be resummed to all orders, there is no 
systematic or unique way to include corrections related to 
the ignored mass of the produced hadron H . In any case, 
"resummed" or "mass corrected" FFs obtained in one 
process should not be used with fixed order expressions 
in other processes. However, it is a key point of global 
QCD analyses of FFs to combine diverse sets of data in 
their determination Only by consistently exploiting 
a large variety of one-hadron inclusive processes one can 
test the underlying assumption of factorization, imply- 
ing that FFs are exactly the same no matter if the seed 
parton is excited from the vacuum like in e + e~ annihila- 
tion or from a nucleon in the case of one-particle- inclusive 
processes in either lepton-nucleon or nucleon-nucleon col- 
lisions. To avoid problems at small z, we follow Ref. [lj 
and simply impose a cut z > z m in — 0.05(0.1) on all 
data with identified pions (kaons) used in our analysis of 
medium modified fragmentation functions. 
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B. Convolution approach for medium modified FFs 

Assuming that hard collinear factorization for one- 
particle inclusive processes is realized also for collisions 
involving nuclei, one can define fragmentation functions 
D^ A (z) specific for nuclear environments characterized 
by their atomic number A. Like standard FFs, the 
medium modified D^j A (z) describe the hadronization of 
a parton i into a hadron H but now in the background 
of a nucleus A. Due to medium-induced final-state soft 
exchanges which happen after the hard partonic scatter- 
ing process, the non-perturbative content of the D*j A (z) 
could differ significantly from the one known for vacuum 
fragmentation functions D^(z). A similar reasoning is 
applied with great phenomenological success in analyses 
of nPDFs that account for medium-induced effects in the 
initial-state p], 03, El • 

Provided that both FFs and nFFs have the same fac- 
torization properties, we can use exactly the same theo- 
retical expressions in computations of measured hadron 
yields for processes involving free or nuclear bounded nu- 
cleons, replacing only the PDFs and FFs by the corre- 
sponding nPDFs and nFFs in the latter case. The global 
analysis of the nFFs, that is the precise determination 
of the medium induced modifications to FFs based on 
the presently available data, proceeds in close analogy to 
those for FFs. The stringent framework of pQCD and 
factorization does not require to make any specific as- 
sumptions on the size and sign of nuclear modifications 
prior to the fit, and the results are entirely determined by 
data. Known nuclear effects on parton densities, as seen, 
e.g., in DIS off nuclei 0, are fully accounted for by choos- 
ing a recent set of nPDFs. The obtained sets of nFFs will 
allow us to test how well the underlying assumption of 
factorization works in practice and can be compared to 
the different ideas and mechanisms for medium modifica- 
tions of quark and gluon FFs proposed in the literature 
0, [1] • We refer the reader to Ref. [lj for detailed expres- 
sions of the relevant cross sections and an outline of the 
numerical fitting procedure. 

Rather than fitting from scratch the nFFs, which 
would take as many parameters as the standard or vac- 
uum FFs, plus several more to represent their nuclear A 
dependence, we choose to relate the Df/ A to the stan- 
dard ones Df at a given initial scale Qo = 1 GeV by a 
convolution approach: 



Df /A {z, Ql) = f ^ W t H (y, A, Q 2 ) Df£ Q 2 ) . 

j z y y 

1 • 1 , r _i • TT7"ff / A . . i • 1 



(4) 



The weight function Wf 1 (y, A, Qq) parameterizes all nu- 
clear modifications and, at the same time, retains the in- 
formation already available on the vacuum FFs for A = 1 . 
Here we take the NLO sets of DSS [lj as reference, which 
provide an excellent global description of hadron yields in 
e + e~, ep, and pp processes. We refrain from performing 
our analysis at leading order accuracy which, at best, can 
give a rough qualitative result. The LO FFs in Ref. [l[ 
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yield a significantly less favorable description of the data. 
Provided that the functional form for Wf 1 (y, A,Qq) is 
flexible enough, one can accommodate the specific de- 
tails of each individual measurement, and the nFFs can 
be extracted as precisely as possible from data. The scale 
dependence of the nFFs in Eq. (j4j is dictated by factor- 
ization and determined by the standard evolution equa- 
tions for vacuum FFs discussed in Sec. Ill Al 

The impact of the weight function W/ 1 (y, A, Qq) in 
Eq. (H]) on the resulting nFFs can be readily understood. 
A simple Dirac delta function 6(1 — y) as weight would 
imply no medium induced effects in the hadronization 
process, while a shift in its argument, i.e., 6(1 — e — y), 
leads to a shift in the momentum fraction z as suggested, 
for instance, by many parton energy loss scenarios 
A more flexible weight function like 

Wf I (y,AQ 2 Q ) = n l y a >(l-yf>, (5) 

can be used to parameterize effects not necessarily related 
to partonic mechanisms, such as hadron or pre-hadron 
attenuation or enhancement, with a great economy of 
parameters. As will be shown below, the A dependence 
of the weight function W^ 1 can easily be included in its 
coefficients, e.g., ni, cti, and Pi in ([5]), by taking them as 
smooth functions of a nuclear property like its volume, 
radius, or mean density. 

We note that convolutional integrals like in Eq. (|3]) are 
the most natural language for parton dynamics at LO and 
beyond, showing up ubiquitously in evolution equations, 
cross sections, etc. They can be most straightforwardly 
handled in Mellin moment space, which is also well suited 
for a numerically efficient global QCD analysis [23[. The 
convolution approach has been demonstrated to be effec- 
tive and phenomenologically successful in the extraction 
of initial-state nuclear effects for PDFs at NLO accuracy 
in Ref. L ll]. For consistency, we take these sets of nPDFs, 
in the following labeled as nDS, in all calculations of cross 
sections relevant for our global analysis of nFFs. 

C. Data sensitive to nFFs 

The gross features of FFs are usually determined from 
very precise e + e~ annihilation data, mainly taken by the 
CERN-LEP experiments. The lack of this source of infor- 
mation for medium modified fragmentation functions is 
evident and significantly complicates their analysis. All 
available probes require a careful deconvolution of nFFs 
from medium effects related to nPDFs, which can be done 
consistently in a global QCD analysis based on factoriza- 
tion. 

One of the most interesting pieces of evidence on 
medium induced effects in the hadronization process 
comes from semi-inclusive deep-inelastic scattering off 
nuclear targets. Such kind of measurements have been 
performed since the late seventies by different collabo- 
rations [5| and in recent years have reached a level of 
precision and sophistication that allows for very detailed 



quantitative analyses. Specifically, the HERMES col- 
laboration has performed a series of measurements on 
deuterium, helium, neon, krypton, and xenon targets, 
with identified charged and neutral pions, kaons, and 
(anti)protons in the final-state 0]. The data are pre- 
sented as distributions in the relevant kinematical vari- 
ables, such as the hadron momentum fraction z and the 
photon virtuality Q 2 , which are used to characterize frag- 
mentation functions, as well as the virtual photon energy 
u, that can be related to the nucleon momentum frac- 
tion x carried by initial-state parton. In addition, data 
are available in terms of p\, the transverse momentum 
squared of the observed hadron. The detailed kinemati- 
cal dependence of the HERMES data 0| , combined with 
the information for different final-state hadrons and tar- 
get nuclei, puts very sharp constraints on the effective, 
medium modified fragmentation functions we wish to de- 
termine. 

In order to minimize initial-state medium induced ef- 
fects, which in a factorized approach should be contained 
in the nPDFs, the data are presented as ratios of hadron 
multiplicities for heavy nuclei A and deuterium (d), 

/ N H (v,Q 2 ,z,p 2 T ) \ 

R H (v O 2 zv 2 )- ± )a (6) 

H A {v,Q ,Z,p T ) - f NH{uQ2 zp 2 TA ■ W 

\ N °{",Q 2 ) ) d 

N H (y, Q 2 , z.p 2 ^) denotes the number of hadrons of type 
H produced in SIDIS, and N e (v,Q 2 ) is the number of 
inclusive leptons in DIS. The cancellation of initial-state 
nuclear effects would be exact in a LO framework if 
medium induced modifications to the PDFs could be 
represented by a single multiplicative factor irrespective 
of the parton flavor. Even though the different parton 
species are known to require different, non-trivial modi- 
fications [llj], these differences get diluted even at NLO 
accuracy. As will be shown below, they cancel in the 
ratio ((6|) to a very good approximation. 

The second crucial piece of evidence on medium in- 
duced effects in the hadronization process comes from 
single-inclusive identified hadron yields obtained in dAu 
collisions at mid-rapidity by the RHIC experiments at 
BNL [Til . ITU [Tr| . These measurements are often seen 
as "control experiments" associated with the program 
of colliding two heavy-ion beams at RHIC to explore 
the properties of nuclear matter under extreme condi- 
tions. However, in the face of the clear evidence of strong 
medium induced effects in the fragmentation pro cess de- 
duced from the SIDIS data @], the dAu data [H [H, [l6| 
also acquire a particular relevance for our analysis. 

Nuclear effects on hadron production in dAu collisions 
are often quantified through comparison to the corre- 
sponding pp spectrum, scaled by the average nuclear 
overlap or geometry function Af . The so defined nuclear 
modification factor R dAu [HI, HE Ell reflects not only 
the medium induced effects on the fragmentation pro- 
cess, but also on the PDFs and includes possible devia- 
tions due to isospin considerations. Af counts the number 
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of underlying inelastic nucleon-nucleon (binary) collisions 
and depends on Glauber model calculations 0, [TH, [l6[ . 
Experimental results for nucleon modification factors and 
invariant cross sections are divided into different central- 
ity classes or presented for the combined minimum bias 
sample. 

Surprisingly, instead of the pronounced hadron atten- 
uation seen in nuclear S1DIS data [0], i.e., R% < 1, the 
nuclear modification factors in dAu collisions show a cu- 
rious pattern enhancement and suppression depending 
on pt [HI, IHl EH ■ The interpretation of these results is 
complicated by the large amount of contributing partonic 
subprocesses with either quark or gluon fragmentation 
into the observed hadron H. Their individual share to 
the hadron yield is strongly correlated with pt as will be 
shown in some detail below. 

Rather than using the invariant cross sections, whose 
theoretical estimates suffer from a sizable dependence on 
the choice of the factorization scale /j, / [24[ , we perform 
the global fit in terms of the ratio of the dAu minimum 
bias cross section to the corresponding invariant hadron 
yield in pp collisions, 



1 Ed 3 a 11 / dp 3 



dA 



2 A Ed 3 a h /dp 3 \ 



(7) 



normalized by the number of participating nucleons. The 
advantage of is that both numerator and denomina- 
tor can be computed consistently within pQCD at NLO 
accuracy using current sets of PDFs [251], nPDFs [TTj | . 
and FFs [l[. In addition, the dependence on \if as well 
as normalization errors associated with inaccuracies in 
determinations of the PDFs and FFs tend to cancel in 
the ratio ([7]). As compared to the related nuclear mod- 
ification factor R^ Au , the ratio is less sensitive to 
estimates based on model dependent calculations. It is 
worthwhile noticing that the DSS set of FFs [l| was ob- 
tained from a global analysis using the same pp collision 
data we use to define the ratios in Eq. . 

The convoluted way in which the information on quark 
and gluon fragmentation is encoded in dAu hadron pro- 
duction data is not an obstacle in a factorized approach 
since we can explicitly compute the relevant cross sec- 
tions in terms of nPDFs and nFFs, as we do for SIDIS. 
In our global analysis, both sets of data will be treated 
simultaneously at NLO accuracy, and the optimum set of 
quark and gluon nFFs will be extracted. In this way we 
can test whether the hadron attenuation found in SIDIS 
off nuclei can be matched with the complicated pattern of 
enhancement and suppression observed in dAu collisions. 

The fact that SIDIS is strongly dominated by quark 
fragmentation suggests that quark nFFs will be sup- 
pressed. Likewise, from the known dominance of gluons 
in pp collisions at RHIC in the relevant low-to-medium 
Pt region [24| . one expects that the main medium in- 
duced effect for them would be an enhancement. In the 
following we describe the details of our global analysis of 
nFFs, present the obtained z and A dependence of quark 



and gluon nFFs for pions and kaons, and show that the 
above, naive expectations are qualitatively correct. 

The obtained nFFs apply to the production of hadrons 
in processes where a nucleus collides with a lepton, a nu- 
cleon, or a very light nucleus like deuterium. For colli- 
sions between two heavy nuclei the distinctive properties 
of the created hot and dense matter presumably may em- 
phasize very different effects in the hadronization which 
perhaps even break factorization. This is beyond the 
scope of this first analysis of nFFs and requires more de- 
tailed studies in the future. 



III. 



RESULTS 



A. Determination of the weight functions 

In order to give a first impression of the main features 
of medium induced effects in the hadronization process 
suggested by the SIDIS @] and dAu [3, [H [l6| data, 
we start with an extremely simple assumption for the 
functional form of the weight functions Wf 1 (y, A, Qq) in 
Eq. ([4]). This will also illustrate the use and convenience 
of the convolution approach outlined in Sec. Ill Bl 

As we have mentioned in Sec.|nCl the HERMES SIDIS 
data show a z dependent hadron attenuation R < I 
that increases with A @. Since the R H do not sig- 
nificantly depend on the different pion charges, H = 
{ir~ , 7T°, 7r + }, it is natural to assume the same medium 
effect on the fragmentation probability for both quarks 
and antiquarks. The observed reduction is most econom- 
ically implemented by a weight function 

Wfiy, A, Ql) = n q S(l -y) + e q 6(1 -e q -y) (8) 

that depends only on two free parameters n q and e q . The 
changes induced in the quark fragmentation functions by 
such a weight are readily interpreted as the superposition 
of two mechanisms. The first term in ([5]) leads to a z 
independent overall reduction of the nFFs relative to the 
vacuum FFs, with n q decreasing from unity with nuclear 
size A. The second term in ([5} takes care of any non- 
trivial z dependence in R H in the most simple way by 
shifting the effective value of z probed in the nuclear 
medium by a small amount e q which grows from zero 
with A. Such a shift is also suggested by models based 
on parton energy loss. 

For gluons, one could try a similar weight function, 
however, low pr pion yields in dAu collisions [l4lll5l 
where gluons are known to dominate in pp processes 
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suggest that the main effect should be an effective en- 
hancement of the fragmentation probability rather than 
a reduction. Therefore n g is most likely to be different 
from n q . Other than in SIDIS, dAu data provide only 
information on the z dependence in a convoluted way, 
i.e., integrated over a certain range in z. Because of this, 
and to a first approximation, we simply take e g — e q . 

The nuclear A dependence of the coefficients n q . g and 
e qi3 parametrizing the naive weight functions W qg in 
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TABLE I: Coefficients parametrizing the naive weight func- 
tions W^ g in Eq. (JSJ) at the input scaie Qo = 1 GeV for dif- 
ferent nuclei A. The three fitted parameters in Eq. ((9} are 
shown in the last line. 



A 


n q 


£ 9 = £ 9 


n g 


He 


0.966 


0.001 


1.015 


Ne 


0.902 


0.002 


1.044 


Kr 


0.745 


0.006 


1.115 


Xe 


0.657 


0.008 


1.155 


Au 


0.550 


0.010 


1.203 




Jn q = -0.0133 


7 E = 0.003 


7„ g = 0.006 



Eq. {8| can be easily implemented by a simple ansatz 

n q = 1 + 7n, ^ 2/3 , 
n 3 = l + 7 „ g A 2 / 3 , (9) 
e q = e g = le A^. 

Such an A dependence could be motivated by the way 
in which the volume of a disk with nuclear radius — 
ro A 1 / 3 , scales with A. In total this results in only three 
free parameters, {7™,, 7n g , 7e}, to be determined by the 
global fit. 

As we will demonstrate in detail below, this very sim- 
ple minded parameterization (in the following labeled as 
nFF*) reproduces to a good approximation the normal- 
ization and general trend of pion yields in SIDIS off nuclei 
and in dAu collisions. For completeness, the obtained co- 
efficients parameterizing the naive weight functions are 
summarized in Tab. |U The ansatz in Eqs. ((8|) and ([9]) 
is, however, not flexible enough to satisfactorily describe 
more detailed features of the data, for instance, regarding 
the x dependence of R\ in SIDIS or the pr dependence 
of the dAu data over the entire range. This is also re- 
flected in the total x 2 (to be defined below) per degree 
of freedom (d.o.f.) of such a fit, which is close to 2. 

These shortcomings suggest that we need to imple- 
ment much more flexible weight functions to fully ex- 
ploit the constraining power of the data. In principle, 
these weights could represent underlying mecha- 

nisms whose consequences would be other than overall 
changes in normalization and shifts in momentum frac- 
tion z, but still preserving factorization. Specifically, we 
adopt the following ansatz for our global analysis of nFFs 
for pions and kaons in Eq. (j4|), 

W q H (y,A,Ql) = n q y a «(l-yf«+n' q 6(l-e q -y), 
W g H (y,A,Q 2 Q ) = ng y"s(l-yfs +n' g 5(l-e g -y), 

(10) 

discriminating between quarks and gluons. With the cur- 
rently available data, no significant improvement of the 
fit is found by introducing an additional weight function 
for antiquarks different from the one used for quarks. In 



addition, the sensitivity of the data on the precise z de- 
pendence of the medium modified gluon fragmentation 
is not sufficient to allow for independent shifts e q and 
e g for quarks and gluons, respectively. Like in our sim- 
ple ansatz in Eq. p|). we set e q — e g , but we allow for 
n' q =/= n ' g . T he parameterizations of the weight functions 
in Eq. (|10|) can be seen as a natural extension of Eq. {SJ, 
where the first term can now have a flexible z depen- 
dence, and free coefficients n! i add extra flexibility to the 
second, energy loss term. 

The nuclear dependence of the coefficients in (TlU)) can 
be again implemented as a smooth function in A, 

£ = \ ( + le A s e, (11) 

where £ = {n q , g ,a q , g , p q>g ,n qg ,e qjg }, and the Ag, 7£ , and 
<5j need to be determined by the fit. Allowing for a com- 
pletely unconstrained A dependence in (fTTjl would require 
too many parameters, which in turn would be only poorly 
constrained by data. Some further guidance is provided 
by the requirement that nuclear effects should vanish as 
A — > 1, and the nFFs should turn into the well-known 
vacuum FFs. Since it is not entirely clear yet that this 
limit is smooth for small A < 4, we do not impose the 
additional constraint that W qg approaches a delta func- 
tion 6(1 — y) as A — ► 1. Nevertheless, in most cases A| 
can be set to zero or unity without damaging the qual- 
ity of the global fit, and the corresponding values of 7^ 
stay very close to zero as required by a smooth limit if 
A — > 1. In addition, some coefficients show a clear pref- 
erence for a nuclear size dependence close to A 2 / 3 . The 
assumption of this simple geometrical behavior with one 
common exponent 6% for all the coefficients in (TTTJ) does 
not change the quality of the fit. In total this leaves 14 
free parameters to be determined by the fit. 

B. NLO analysis of pion nFFs 

The free parameters in Eqs. (TTU|) and (TTTj) , as well as 
for the simplified fit presented in the previous Subsection, 
are determined by a standard \ 2 minimization for N data 
points, where 

i— 1 1 

Ei is the measured value of a given observable, 6Ei the 
error associated with this measurement, and Tj denotes 
the corresponding theoretical estimate for a given set of 
parameters. As often, we take statistical and systemat- 
ical errors in quadrature in 5Ei. Since this is the first 
attempt of an extraction of nFFs from a global QCD 
analysis, we refrain from a more sophisticated treatment 
of experimental uncertainties. In any case, full informa- 
tion on error correlation matrices is not available for most 
of the data sets being analyzed here. 

The results of a fit based on the ansatz for the weight 
functions W q g presented in Eqs. (TTUj) and (jlip is sum- 
marized in Tabs. |TT] and IIIII In Table [IT] we collect the 
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TABLE II: Coefficients parametrizing the weight functions W q>g in Eq. (|10[) at the input scale Qo = 1 GeV for different nuclei 
A. The 14 fitted parameters A^, 75, and 8$ in Eq. (|11[) are given in the bottom half of the table. 



A 


n' q 


£ 9 = £ 9 


n q 


a q 


P, 


n' g 


n g 


a 9 




He 


0.949 


0.002 


0.024 


24.18 


26.80 


1.042 


-0.131 


19.50 


46.24 


Ne 


0.863 


0.006 


0.065 


23.54 


26.13 


1.114 


-0.351 


19.21 


41.78 


Kr 


0.668 


0.015 


0.157 


22.10 


24.62 


1.276 


-0.849 


18.56 


31.72 


Xe 


0.564 


0.020 


0.206 


21.33 


23.81 


1.362 


-1.115 


18.21 


26.33 


Au 


0.439 


0.026 


0.265 


20.41 


22.85 


1.466 


-1.433 


17.79 


19.90 




1 








24.56 


27.20 


1 





19.67 


48.88 


75 


-0.022 


0.001 


0.010 


-0.161 


-0.169 


0.018 


-0.056 


-0.073 


-1.126 


S( 


0.615 


0.615 


0.615 


0.615 


0.615 


0.615 


0.615 


0.615 


0.615 



TABLE III: Data sets included in the NLO global analysis of 
pion nFFs, the individual x 2 values for each set, and the total 
X 2 of the fit. 



Data Data 



Experiment 


A 


H 


type 


points 


x 2 


HERMES [6] 


He,Ne,Kr,Xe 


7T + 


z 


36 


39.3 






7T~ 


z 


36 


23.0 






7T° 


z 


36 


27.4 






7T+ 


X 


36 


69.4 






7T~ 


X 


36 


55.4 






7T° 


X 


36 


49.7 






7T+ 


Q 2 


32 


21.0 






7T~ 


Q 2 


32 


27.1 






7T° 


Q 2 


32 


34.7 


PHENIX [14] 


Au 


7T° 


pT 


22 


13.7 


STAR (prel.) [16] 


Au 


7T° 


[>T 


13 


12.8 


STAR [15] 


Au 


7r± 


Pt 


34 


22.5 


Total 








381 


396.0 



values of the coefficients parametrizing the weight func- 
tions Wg g in Eq. (fTT))) for different nuclei and the fitted 
parameters A^, 7^, and 8% in Eq. (fTTjl . Tablc lllTl shows the 
partial contributions to x 2 f° r each set of data included 
in the fit. We obtain an overall x 2 — 396.0 for 381 data 
points included in the analysis with 14 free parameters, 
resulting in an excellent x 2 /d.o.f = 1.08. Recall that the 
parameters of the simplified, three parameter fit based 
on Eqs. © and ©, which yields x 2 /d.o.f . ~ 2, were 
already presented in Tab. Q] 

The obtained smooth A dependence of the coefficients 
parametrizing the weight functions W^ g in Eq. (fT0|) is 
illustrated in Fig. [TJ Although not explicitly enforced in 
Eq. (jlip . n' g approaches unity and both n g>g and e' q g 
tend to zero as A — ■> 1, as required by the vanishing 
of nuclear effects in that limit. As expected from the 
qualitative discussion of the data in Sec. Iff CI the pat- 
tern of medium induced modification is rather different 
for quarks and for gluons. Most importantly, the nor- 



malizations n' q and n' q of the Dirac delta term in the 
convolution weights (|10[) have opposite signs, leading to 
suppression for the quark and enhancement for the gluon 
nFFs with respect to the vacuum FFs. The first term on 
the right-hand-side of Eq. (jTUJ) also shows an opposite 
trend for quarks and gluons but influences mainly the 
small z behavior of the resulting D q , A and D*, A which 
will be presented at the end of this Subsection. 

In Fig. [2j the measured multiplicity ratios R\ for 
charged and neutral pions Q are shown as a function 
of z, x, and Q 2 for different target nuclei A. To demon- 
strate the significance of medium modifications in the 
hadronization process, the dashed lines correspond to cal- 
culations of R A at NLO accuracy using nPDFs from nDS 
p"l1 | but standard DSS fragmentation functions [l[. All 
computations of multiplicities on a deuterium target in 
the denominator of Eq. ^ are performed with MRST 
PDFs [H and DSS [l| FFs. As anticipated in Sec. [EC] 
initial-state nuclear effects approximately cancel in i?^, 
and the results computed in this way are very close to 
unity for all kinematic distributions, in sharp contrast to 
data. Very similar results are obtained if other current 
sets of nPDFs [H, d are used. 

Notice that despite cancellations of the initial-state ef- 
fects associated to the nPDFs, the measured medium in- 
duced modifications of the pion multiplicities can be as 
large as a fifty percent effect for the heavier nuclei. They 
show a non-trivial z and x dependence, increase with the 
nuclear mass number A, and are most conspicuous for 
larger hadron energy fractions z and larger x, i.e., smaller 
v. The dependence on Q 2 , displayed in the bottom panel 
of Fig. is comparatively flat but noticeable. Several 
models proposed to estimate medium induced effects in 
the hadronization process reproduce some of the features 
of the data. However, the full kincmatical dependence of 
R A (u, Q 2 , z,Pt) is still a challenging issue; see Refs. 
and references therein. 

The results of the global analysis of nFFs based on the 
convolution approach are shown as dotted and solid lines 
for the naive three parameter and the refined ansatz for 
the weight functions W qg introduced in Sec. IIII Al re- 
spectively. Both fits give a much superior description 
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FIG. 1: A dependence of the coefficients parametrizing the weight functions W£ g in Eqs. (fTO]) and (fTTj) . 



of the full kinematic dependence of the HERMES data 
than an approach which ignores final-state nuclear ef- 
fects. Even the simple ansatz is doing surprisingly well 
in reproducing the general trend of the data, with ten- 
sions mainly for the x and z differential yields and larger 
nuclei. It is interesting to notice that there seems to be 
no visible conflict between the standard Q 2 dependence 
assumed for the nFFs in our fit and the data. In this re- 
spect, there have been many interesting suggestions and 
model dependent calculations at the LO level, motivat- 
ing the use of medium modified evolution equations, see 
[29 ] and references therein. However, in the range of Q 2 
covered by present SIDIS data, there is no evidence for 
a significant departure from standard timelike evolution 
equations and kernels in Eq. |T]) and ([3]), respectively. 

Even though the full NLO framework for the px de- 
pendent hadron yields in SIDIS processes is available 
[271 ] , and very precise experimental studies have been pre- 
sented recently we choose to work only with px in- 
tegrated R*[ for the time being. The dependence of the 
data on this variable is rather weak, and the px values 
accessible so far are at the limit of a perturbative treat- 
ment. 

In Figs. G2 [U and [5] we show both the invariant dAu 
cross sections for neutral and charged pion production at 
RHIC as a function of px, and, in the lower panels, the 
ratios R„ to the pp yields, defined in Eq. ((7]). Due to 
the steep fall of the cross sections with px over several 
orders of magnitude in the range shown in Figs. 00 it 



is hard to see any differences between the data and the 
NLO calculations based on standard or medium modified 
FFs. They become clearly visible, however, in terms of 
the ratios R% and can be as large as about 20%. 

At variance with R\ in SIDIS, where initial-state 
medium effects, and hence the differences associated to 
the choice of nPDFs, cancel to a large extent, they may 
remain significant for dAu data. Thus, estimates of the 
medium induced modifications for fragmentation func- 
tions could in principle depend on the choice of nPDFs. 
To estimate their impact on the dAu cross sections and 
on the ratios i?^, we use two different sets of nPDFs, nDS 
fPH and EPS [lj], along with the well-known standard 
FFs of DSS for our calculations shown in Figs. [Mil As 
can be seen, both results do not differ too much and fall 
short of describing the data satisfactorily, leaving room 
for improvement due to medium induced effects in the 
hadronization process. The EPS set reproduces the trend 
of the data better than nDS, in particular for neutral pi- 
ons shown in Figs. [3] and [4] which extend to somewhat 
larger values of px than the charged pion data in Fig. [5] 
Notice, however, that the dAu data were included in the 
EPS analysis of nPDFs, but assuming nuclear effects in 
the hadronization process to be negligible [13]. Although 
the modification introduced in the EPS nPDFs helps to 
describe the dAu data better [13], neglecting final-state 
nuclear effects is clearly not advisable in view of their sig- 
nificant impact on R\ in SIDIS demonstrated in Fig. [2j 

The results of our global fit of nFFs are shown in 
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FIG. 2: JJ^in SIDIS for different nuclei in bins of z (upper panel), x (middle panel), and Q 2 (bottom panel) as measured by 
HERMES [a|. The solid lines correspond to the results of our optimum fit for nFFs using the nDS medium modified parton 
densities The corresponding fit based on the simple nFF* ansatz in Eq. |(5J is shown as dotted lines. The dashed lines are 
estimates assuming the nDS medium modified PDFs but standard DSS vacuum FFs [J]. 



Figs. [3][5] as solid and dotted lines, corresponding to our 
optimum and simplified ansatz for the weight function 
W£ g introduced in Sec. IIII A[ respectively. Here, the 
naive three parameter ansatz for W£ g fails to reproduce 
the pt dependence of the dAu data, and the greater flex- 
ibility of Eqs. (fT0|) and (lll|) is clearly needed and leads 
to a significant improvement of the fit. The simultaneous 
description of SIDIS and dAu data requires to have the 
correct balance between quark and gluon contributions 
in the fragmentation process, which is strongly and, 
hence, z dependent. 



An important difference between SIDIS and dAu data 
is that in the latter case the cross sections sample con- 
tributions from a wide range in z. Consequently, the 
deconvolution of the medium induced effects is less trans- 
parent. In order to provide a better insight into the sen- 
sitivity of the RHIC measurements to the fragmentation 
process, we show in Fig. [5] (a) the mean value of the 
hadron's fractional momentum z probed in pp and dAu 
collisions as a function of pr- There are several ways to 
estimate an average (z). We define it in the standard 
way by evaluating the convolutions in the factorized ex- 
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P T [GeV] 

FIG. 3: Upper panel: comparison of the PHENIX data for 
neutral pion production in dAu collisions at mid-rapidity [l4| 
with NLO estimates obtained with various combinations of 
nPDFs, FFs, and nFFs. The solid and dotted line correspond 
to the results of our optimum and simple three parameter fit 
for nFFs, respectively, usin s the nDS nPDFs hj. Results 
based on standard DSS FFs [1| are shown as dashed and dot- 
dashed lines for nDS [U and EPS [H nPDFs, respectively. 
Lower panel: same as above but now for the ratio RZ defined 
in Eq. ©. 



pression for the pt dependent cross section [24j with an 
extra factor of z in the integrand, divided by the cross 
section itself [28j . i.e., schematically we use 



(z) 



j dz ; 



da" 
dzdpi 



J dz dzdpT 



(13) 



Here, da H / dzdpx contains the appropriate convolutions 
of the parton densities and fragmentation functions with 
the partonic hard scattering cross sections. 

Panel (b) of Fig. [5] shows the individual (z) for quark 
and gluon fragmentation processes, simply referring to 
the contributions in the dAu cross section proportional to 
either D^ Au or D^ Au . Beyond the LO, this separation 
involves some arbitrariness and depends on the choice of 
the factorization scheme. In addition, primary partons 
created in the hard scattering may radiate off secondary 
partons of a different species which in turn fragment into 
the observed hadron. Figure [6] (c) shows histograms of 
the z distribution for three representative values ofpx- In 
panel (d), we present the relative contributions of quark 
and gluon fragmentation processes to the ir° production 
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FIG. 4: The same as in Fig. [3] but now for neutral pion data 
obtained by STAR 



cross section in pp and dAu collisions. 

As can be seen in Fig. El RHIC pp and dAu data 
are mainly sensitive to fairly large values of the momen- 
tum fraction taken by the hadron H , with (z) slightly 
increasing with px- However, as panel (c) shows, the 
cross section samples contributions over a broad range in 
z, starting at about z ~ 0.2. Notice that below about 
Pt = 1-5 GeV, the tail in the z distribution becomes sen- 
sitive to values z < 0.1, where the concept of fragmen- 
tation functions breaks down due to finite hadron mass 
effects and the singular behavior of the timelike evolution 
kernels. It is also worthwhile mentioning that the values 
of z to which the cross sections are most sensitive to, i.e., 
(z), depend, of course, on the actual shape of the FFs 
and nFFs assumed in the analysis of pp and dAu collision 
data, respectively. Since we anticipate sizable differences 
between them, the ratios defined in Eq. ([7]) actually 
sample the nuclear and the vacuum fragmentation func- 
tions at slightly different values of z, cf. Fig. [5] (a). This 
can have quite some effect on the ratios R£ in regions 
where the fragmentation functions vary rapidly with z. 

Fig. [6] (d) demonstrates that for pp collisions gluon 
fragmentation is clearly the dominant production mech- 
anism at low values of pt- Quark fragmentation con- 
tributions increase with pt, and cross the level of 50% 
at pt — 10 GeV. Again, the relative balance between 
quark and gluon contributions in dAu collisions will de- 
pend on the extracted medium induced modifications. 
In our analysis, pion production is dominantly driven by 



11 




P T [GeV] p T [GeV] 



FIG. 5: The same as in Figs. and U but now for charged pion data from STAR [l5l ]. 




FIG. 6: (a) The mean value of z probed in pp and dAu collisions as a function of pr- (b) The mean value of z in dAu collisions 
for both quark and gluon contributions, (c) Histograms of the z distribution for three representative px values, (d) The relative 
contributions of quark and gluon fragmentation processes to the 7r° production cross section in pp and dAu collisions. 
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FIG. 7: The resulting medium modified NLO fragmentation functions for quarks and gluons into neutral pions at Q 2 — 10 GeV 2 
for various different nuclei (right panels). The left panels show the corresponding ratios R~ „ to the standard DSS FFs as defined 
in Eq. 



gluon fragmentation up to the highest values of pt, about 
15 GeV, currently accessible in experiment. As we shall 
see below, this is due to the nuclear suppression of D^, A 

on the one hand, and enhancement of D H , A on the other 
hand. 

The resulting modified fragmentation functions and 
the corresponding ratios to the standard DSS FFs, 

H " - D»(z,Q 2 ) ' Ua ~ Df(z,Q*) ' [ > 

can be found in Fig. [7] at Q 2 — 10 GeV 2 for various nu- 
clei. As expected from the qualitative discussion of the 
data in Sec. Ill C[ the pattern of medium induced modifi- 
cations is rather different for quarks and for gluons. The 
dominant role of quark fragmentation in SIDIS leads to a 
suppression, i.e., < 1, increasing with nuclear size A 
as dictated by the pattern of hadron attenuation found 
experimentally, see Fig. O The enhancement of hadrons 
observed in dAu collisions for pt < 10 GeV, see Figs.[3][5l 
along with the dominant role of gluon fragmentation at 
low values of pr explains that R* > 1 for z > 0.2. Below 
z ~ 0.2, where all the data we analyze have very little 
or no constraining power, both quark and gluon nFFs 
drop rapidly. For the time being, the behavior in this re- 
gion could easily be an artifact of the currently assumed 
functional form for the weig hts W£ g in Eq. (TU|). 



We note that for z > 0.4 the bulk of the medium in- 
duced effects for both quarks and gluons can be accom- 
modated by the second term in W^ g proportional to a 
Dirac delta function, which normalization coefficients n' q 
and n' g have to have opposite signs, reflecting suppres- 
sion and enhancement, respectively. The other term in 
Wq g introduces only a small correction in this region, 
but becomes dominant if z < 0.4, especially for D^, A . 

We wish to stress, that even with the limited amount 
of data available at present, the successfully performed 
global analysis of nFFs provides a first non-trivial indi- 
cation that the assumed factorization of long- and short- 
distance physics works amazingly well also for SIDIS off 
nuclei and dAu collisions at RHIC. The found suppres- 
sion R* < 1, as imposed by SIDIS data, is fully compat- 
ible with the complicated pattern of enhancement and 
attenuation of hadron yields in dAu collisions at RHIC. 
The non-negligible role of quark fragmentation in dAu 
collisions at moderate and large values of px, along with 
the larger nuclear effects for D^ A than for D^ A at large 
values of z, explains that R£ < 1 for pt > 10 GeV. In 
addition, even though SIDIS data are dominantly sen- 
sitive to quark fragmentation, the observed non-trivial 
Q and x dependence of the SIDIS multiplicity ratios 
shown in Fig. [21 depends on the medium induced mod- 
ifications of the gluon fragmentation function, which in 
turn is mainly constrained by dAu data. The correlation 



13 




0.5 



+ 



10 



P T [GeV] 



FIG. 8: As in Fig.[3]but now separated into different centrality 
classes ranging from the most central one, 0—20%, to the most 
peripheral one, 60 — 88%. The theoretical estimates obtained 
for minimum bias events are rescaled by the ratio C denned 
in Eq. 1|15[1 to account for the different centrality classes. 
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FIG. 9: As in the lower panel of Fig.|8]but now in terms of the 
experimentally determined nuclear modification factors RdAu 



of the x dependence of SIDIS multiplicity rates with the 
z dependence of the fragmentation functions is induced 
by the NLO hard scattering coefficient functions which 
depend in a non-trivial way on both x and z; see, e.g., 
Ref. [H. 



minimum bias (m.b.) sample, i.e., 



C = 



Ed 3 a*/dp 3 
Ed 3 a" / dp 3 



dAu \ 
m.b. I 
dAu ' c.c. 



(15) 



C. Dependence on different centrality classes 

All the data on hadron yields in dAu collisions we have 
discussed and analyzed so far correspond to minimum 
bias events. Experiments also present their results di- 
vided into different centrality classes depending on how 
central or how peripheral the collision is in impact pa- 
rameter space. In that respect, the nFFs obtained in our 
global analysis correspond to some sort of average of the 
medium induced effects seen for the different degrees of 
geometrical overlap between the deuteron and the gold 
nuclei. One could expect that this average underesti- 
mates the medium induced effects when the collisions 
are more central and overestimates them for the more 
peripheral ones. 

To estimate the possible impact of different centrality 
classes on the extraction of nFFs, we scale our results 
obtained for the minimum bias cross sections shown in 
Figs. [3][S1 by the average ratio between the measured 
cross section for a given centrality class (c.c.) and the 



The ratio (fT5"]) is a simple estimate of the fraction of 
events selected by a given centrality cut. 

As a representative example, we show in Fig. [5] the 
comparison between the PHENIX data for four different 
centrality classes [lj| and the corresponding minimum 
bias cross section computed with our nFFs and rescaled 
by C. As expected, our nFFs slightly overestimate the 
data for the most peripheral events denoted by 40 — 60% 
and 60 — 88%. The agreement with the more central 
events is, however, very good, and we do not find the ex- 
pected underestimate of nuclear effects. As for the mini- 
mum bias cross sections shown in Figs.[5][5l estimates ob- 
tained with standard vacuum FFs do not reproduce the 
trend of the data very well, except for the most periph- 
eral events. Results for pion yields in different centrality 
classes obtained by the STAR experiment are very simi- 
lar to those presented in Fig. [8] and hence not shown. For 
completeness, Fig. Ogives a comparison similar to those 
presented in the lower panel of Fig. [5] but now in terms 
of the nuclear modification factors. 
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TABLE IV: Coefficients parametrizing the weight functions W^ g in Eq. (|10|> at the input scale Qo = 1 GeV for different nuclei 
A. The 13 fitted parameters A^ and 7^ are given in the bottom half of the table. Note that the are chosen as in Tab. [TTJ 



A 


n q 


c 


n q 


Ctq 


Pi 


n„ 


n 9 




Pa 


He 


0.957 


0.003 


0.024 


27.95 


30.16 


1.038 


-0.107 


17.75 


39.37 


Ne 


0.884 


0.007 


0.065 


27.37 


29.55 


1.103 


-0.288 


17.40 


35.33 


Kr 


0.719 


0.017 


0.156 


26.08 


28.17 


1.248 


-0.695 


16.60 


26.22 


Ze 


0.630 


0.022 


0.205 


25.38 


27.43 


1.326 


-0.913 


16.17 


21.34 


Au 


0.525 


0.028 


0.264 


24.55 


26.55 


1.419 


-1.174 


15.67 


15.52 




1 








28.29 


30.52 


1 





17.96 


41.76 


7€ 


-0.0185 


0.0011 


0.0102 


-0.1451 


-0.1543 


0.0163 


-0.0456 


-0.089 


-1.020 



D. NLO analysis of kaon nFFs 

Even though the available experimental information 
on kaon production in a nuclear environment is much 
more limited than in the case of pions, it is interesting to 
study to what extent the pattern of nuclear modifications 
found for pions is similar to those for kaons. It should be 
noted that the standard kaon FFs, which will provide the 
baseline to analyze medium induced effects, still suffer 
from sizable uncertainties 

In a scenario where the medium induced modifications 
to the hadronization process are dominated by partonic 
mechanisms, and to a first approximation in a convo- 
lutional approach, one expects a very similar pattern of 
medium modifications for pions and for kaons. If interac- 
tions of the produced hadron or intermediate pre-hadrons 
with the nuclear environment are important, medium in- 
duced effects in the production of pions and kaons can 
be significantly different. 

Due to the relative scarcity of data with identified 
kaons, especially in the case of dAu collisions where most 
of the data are taken at pr values below the reach of 
perturbative QCD methods, we cannot proceed with our 
global analysis as we did pions. The coefficients of a com- 
pletely flexible parametrization for the weights W^ g like 
in Eq. (p~0|) would not be well constrained by the data. 
Therefore, our strategy is as follows: we start by impos- 
ing a more constrained ansatz for the weights W^g where 
we assume some of the features of the medium induced ef- 
fects found for the pion nFFs. More specifically, we keep 
the same functional form for the weights as in Eq. (TIT))) , 
with the same nuclear A dependence as in Eq. (fTTj) . We 
also take over the simplifying assumption e q = e g , which 
we have chosen for pion nFFs in Sees. IIII Al and IIII Bl 
after checking that more flexible options did not lead to 
a significant improvement of the quality of the fit. In 
addition, we try to set the exponents 5^ governing the 
A dependence in Eq. (jTTJ) to the values preferred by the 
pion data. 

Next, we can assess how much variation of the other 
parameters relative to the results for the pion nFFs listed 
in Tab. |TT] is required to reproduce the main features of 
the kaon data. We find that allowing for an up to 20% 



TABLE V: Data sets included in the NLO global analysis of 
kaon nFFs, the individual x 2 values for each set, and the total 
X 2 of the fit. 

Data Data 



Experiment 


A 


H 


type 


points 


x 2 


HERMES [6] 


He,Ne,Kr,Xe 


K+ 


z 


36 


61.3 






R- 


z 


36 


91.7 






K+ 


X 


36 


67.2 






K~ 


X 


36 


128.5 






K+ 


Q 2 


32 


38.4 






K- 


Q 2 


32 


50.0 


STAR [15] 


Au 


K+ 


PT 


5 


5.4 




Au 


K~ 


Pt 


5 


4.6 


Total 








218 


447.1 



variation yields a more than reasonable agreement with 
all available data sets, with most of the parameters stay- 
ing within a 10% variation of their counterparts describ- 
ing the pion nFFs. In Tab. IIVI we list the values of the 
coefficients parametrizing the weight functions W^g in 
Eq. (flTJ|) for different nuclei and the 13 fitted parameters 
and 75 in Eq. (JTTJ) ; the are fixed to the values given 
in Tab. \H\ 

The overall quality of the fit is summarized in Tab. [V] 
The total x 2 of the fit is 447.1 for 218 data points in- 
cluded in the global analysis, resulting in \ 2 /d.o.f. = 2.2. 
Also given are the individual contributions to \ 2 for each 
set of data included in the fit. Even though the qual- 
ity of the fit is not as good as in the case of pions, it 
reproduces the data well within their uncertainties and 
suggests a close similarity between the medium induced 
modifications for pion and kaon fragmentation functions. 
As it could be expected, the most significant variations 
are found in the parameters related to the gluon fragmen- 
tation Dg, A , which are only poorly constrained by the 
scarce dAu data. Note that the largest individual contri- 
bution to x 2 stems from the SIDIS data with identified 
K~ where already the vacuum FFs show some tension 
with the experimental results [l|. 

In Fig. [TD] we compare the result of our fit with the 
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FIG. 10: As in Fig. [2] but now for identified K d 



measured SIDIS multiplicity ratios R~j[ for charged kaons 
as a function of z, x, and Q 2 for different target nuclei 
A. As can be seen, the agreement with the K + data 
is reasonably good with some problems in reproducing 
the x distributions for heavier nuclei. Discrepancies are 
somewhat larger for K~ data as we have already men- 
tioned. For comparison, we show again the results of 
a theoretical calculation based on standard vacuum FFs 
(dashed lines). As for the pion multiplicity ratios shown 
in Fig. [21 not even the trend of the data can be repro- 
duced by ignoring medium induced modifications to the 
hadronization process. Finally, Fig. QT] shows both the 
invariant dAu cross section for charged kaon production 
as a function of pt, and, in the lower panels, the ratios 
Re to the corresponding pp yields. One should notice 
the limited range in pt covered by the presently available 
data which is certainly at the borderline where pQCD is 
applicable. Nevertheless, our fit describes the data well 
except for the lowest bin in pr. 

The availability of both more precise kaon data and 
more accurate sets of vacuum FFs in the future, will help 



to further investigate the interesting and non-trivial close 
relation between the nuclear modifications for pion and 
kaon fragmentation functions found in our analysis. 

IV. SUMMARY AND CONCLUSIONS 

We have investigated the possibility of using a fully 
factorized approach at NLO accuracy of pQCD, similar to 
those established for analyses of standard FFs and PDFs, 
to describe medium induced effects in the hadronization 
process. 

To this end, we have explored the concept of modified 
fragmentation functions which effectively accounts for all 
the distortions in the production of hadrons caused by the 
nuclear environment. These novel distributions, which 
we denote as nFFs, are completely analogous to nuclear 
PDFs, and, combined with them, allow one to treat a 
large class of hard reactions where a nucleus collides with 
a lepton, a nucleon, or a very light nucleus like a deuteron 
in a consistent pQCD framework. 
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We have performed combined NLO fits to the avail- 
able pion and kaon production data in SIDIS and in dAu 
collisions to reveal the non-perturbative features of the 
nFFs, exploiting the well established potential of pQCD 
to describe hard scattering processes and the wealth of 
information compiled in PDFs and nPDFs over the past 
years. We choose to relate the nFFs to the well-known 
standard FFs in a convolutional approach with a great 
economy of free fit parameters, including a smooth de- 
pendence on the nuclear size A. 

The obtained parameterizations of nFFs for pions and 
kaons [30j rather accurately reproduce the data and their 
nuclear A dependence. Our results provide support to 
the idea that conventional factorization of short and long- 
distance physics effects works to a good approximation 
also in the nuclear environments studied in our analyses. 
The found pattern of medium induced modification is 
rather different for fragmenting quarks and gluons, where 
we find suppression and enhancement, respectively, com- 
pared to the vacuum fragmentation functions. Since our 
nFFs parametrize the available data without resorting 
to a certain model or underlying mechanism for the ob- 



served medium modifications, they may help to further 
our understanding on hadronization in a nuclear envi- 
ronment. Predictions based on the observed pattern for 
quark and gluon nFFs can be tested by upcoming data 
from BNL-RHIC and CERN-LHC and perhaps in the 
future at an electron-ion collider like the EIC project. 
This will help to understand the limits of the factorized 
approach which is expected to be only an approxima- 
tion and may receive corrections beyond the leading twist 
level or beyond a pQCD approach. 
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